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" Abstract 

9+ 

The time evolution of low energy spin states of a single molecular magnet in a local electric field 
is investigated. The decoherence of the driven single molecular magnet weakly coupled to a thermal 
bosonic environment is analyzed by the second-order time-convolutionless non-Markovian master 
equation. If the characteristic time of the system is much smaller than the correlation time of the 
environment, the analytical expression of the reduced density matrix of the system is obtained. 
The non-Markovian dynamics of the spin states at low temperatures is induced by the memory 
effects in the decay rates. The non-Markovian oscillation of the Bloch vector is presented. The 
quantum decoherence can be effectively restrained with the help of the reasonable manipulation of 
the environment spectral density function and local electric field. The influence of the dissipation 
on the pointer states are investigated by the von Neumann entropy. The pointer states can be 
selected by the environment. 
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I. INTRODUCTION 



The unavoidable interactions of all open quantum systems with environments often result 



in the dissipation and decoherence l[ [2| . Due to the exchange of energy and information 
between the system and the environment, the non-Markovian dynamics of quantum states 



always occur in the realistic experimental systems 



5). Recently, much attentions have been 



paid to the control of the decoherence of many-body quantum systems such as spin 

clusters and single molecular magnets js]. As a class of systems with rich quantum properties, 
single molecular ma gne ts at low energies can serve as a large-spin system or a collection of 
interacting spins j9-ll|. These solid quantum spin systems are considered as promising 



carriers of quantum information 



12|, Ll3[] - Single molecular magnets with antiferromagnetic 



spin couplings can provide low-energy states for performing quantum logic gates 



Ij4l6|. 



Quantum decoherence is manifested when single molecular magnets are coupled to a spin 



bath or thermal reservoir 



flfl 



17j-2(}. The dissipation and decoherence always depend 



on the properties of the environment which can be described by a certain spectral density 



function 



21 



Lri24j. Therefore, a reasonable quantum manipulation method is necessary. At 

present, chemical manipulation can offer an efficient way to engineer intermolecular couplings 

I I 

and allow for interactions between qubits [251 ] . The decoherence from the chemical control 



cannot be easily eliminated because of the permanent interactions with the surrounding 



171 ]. Simultaneously, the most straightforward and conventional way is to adopt an external 
magnetic field produced by electron spin resonance pulses Although the decoherence 
of single molecular magnets can be suppressed by strong magnetic fields, it is preferable to 
apply electric fields that are controllable and suitable on very small spatial and temporal 
scales. It is possible to apply time-dependent str ong electric field close to single molecular 



magnet via a scanning tunnel microscopic tip 26|, |27j. It has been found out that an electric 



field can be coupled to low-energy spin states of different chirality due to the absence of 
spin inversion symmetry in some single molecular magnets, such as C113 [26] , V15 28], 
Co 3 Q, Dy 3 3 and Mn l2 QQ, etc. The effective spin electric coupling relies on the 
detailed structure of single molecular magnets at low energies 3l|. Moreover, the use of 
microwave cavities can contribute to the indirect generation of fully controllable and long- 
range interaction between any two molecular magnets. This scheme based on electric-field 
local control can open up the possibility of scalable solid quantum information processing. 
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Thus, the control of the decoherence of single molecular magnets driven by local electric 
field in low-temperature environments needs to be further studied. 

This paper is organized as follows. In Sec. II, the time evolution of low-energy spin states 
in an electrically driven single molecular magnet without the spin inversion symmetry is 
studied. Using the general low-energy spin Hamiltonian with effective spin-electric coupling, 
the time-convolutionless non-Markovian master equation is derived under the assumption 
of weak interactions with a thermal bosonic environment. If the characteristic time scale 
of the system is much shorter than the correlation time of the environment, the reduced 
density matrix of the low-energy spin states can be obtained analytically. As an example, 
the effects of the low temperatures on the decaying rates of Cu% are calculated. In Sec. Ill, 
the non-Markovian decoherence of spin states is described by means of the Bloch vector. 
The traditional Lorentzian spectral density function of the environment is considered. The 
selection of the pointer state measured by the von Neumann entropy is also investigated. 
Finally, a simple discussion concludes the paper. 



II. MODEL OF DECOHERENCE IN THERMAL BOSONIC ENVIRONMENT 



Molecular magnets have clear features of coherent behavior and a variety of effective low- 
energy spin Hamiltonian is used for encoding qubits and implementing spin-based quantum 
computation. A local electric field e{t) can be coupled to low-energy spin states of opposite 
chirality because of the effective spin electric interactions 3l|. It is shown that both spin- 
orbit interactions and the absence of spin inversion symmetry can induce the the electric 
dipole matrix element d which is an important quantity in the effective spin-electric coupling. 
The strength of spin-electric coupling can be calculated by means of ab initio methods. For 
a single molecular magnet with an effective spin-electric coupling, the effective low-energy 
spin Hamiltonian is written as 

Heff — Hq + H e , (1) 

where H is the low-energy spin Hamiltonian without the electric field and the effective 
spin-electric coupling is given by H e = d ■ e = Y2i ' Here e is the electron charge and r*j 
denotes the coordinate of the ith electron in the spin structure of single molecular magnet. 
In the following discussion, a simple case of triangular spin-| molecular magnet Cu% 26] is 
considered. In the low-energy Hilbert subspace spanned by two opposite chiral spin states 
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{\C = ±1,M = i)}, the effective spin Hamiltonian in the presence of effective spin-orbit 
interaction oj so and planar electric field e{t) can be expressed as 

H eff = l -u so C z + d ■ e{e-^ t+ ^C + + S^C-]. (2) 

The chirality operator is C z = t^^i • (s*2 x s* 3 ) where s*j denotes the spin operator of zth site. 
The operators C± reversing the chirality of the spin states satisfy that C±\C = =)=1,M) = 
\C = ±1, M) and C±\C = ±1, M) = 0. The chirality operator behaves similar to the spin 
operator in the chiral space. Here \C = ±1, M) are the simultaneous eigenvectors of C z and 
total spin operator S z = sf with the corresponding eigenvalue C and M. The ground 
state M — | is taken into account. The parameter u> represents the frequency of the electric 
field, d = \d\ describes the strength of the electric dipole. With the variation of the initial 
angle between the field e and the vector r*i2 from site 1 to site 2, we can reasonably adjust 
the angle a = 



ail. 

A thermal bath can be described by an infinite chain of quantum harmonic oscillators. 
In the rotating frame with the electric field frequency u, the total Hamiltonian of the open 
system coupled to the environment can be written as 

H = H e ff + He + Hj. (3) 

The effective spin Hamiltonian of the molecular magnet is transformed to 

H eff = ^(A so C z + d-eC x ) J (4) 

where A so = u so — ui and the chirality operator C x = \{C + + C_). The Hamiltonian of the 
thermal environment is given by 

H E = Y J ^b% (5) 

j 

where bj and fot are the annihilation and creation operator in the Hilbert space of bosonic 
environment. The last term in Eq. (3) denotes the weak interaction between the system 
and the environment, and can be written as 

Hi = E^e-^K + g*e^%C + ). (6) 
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Here the rotating wave approximation is adopted. The weak coupling \gj\ = d ■ €j where ej 
represents the magnitude of the electromagnetic field of the jth mode with the frequency Uj. 



/ dhTTEiHft), [H'jih), pit) ® PE }}, (7) 
Jo 



For the convenience, the Planck constant h is set to be one. In the interaction representation, 
the decoherence of the low-energy spin state p(t) can be approximately given by the second- 
order time-convolutionless master equation, 

dp{t) 
dt ~ 

where H'j(t) = e ^#e//+tfs)# je -i*(# e //+#s). The notation Tr^ is the partial trace over the 
freedom of the environment. It is assumed that the initial product state of the total sys- 
tem is ptoi(O) = jo(0) <S> Pe where pE = exp(— He / krT) /Tr\exp(— Hr/kbT)] is the thermal 
equilibrium state of the environment and satisfies that Tr[Hj(t)pE} = [24]. The low tem- 
perature condition of kbT < u> so is considered here. For the convenience, the Boltzmann 
constant k b is also assumed to be one. 

To simplify the analytical calculation, the effective spin Hamiltonian can be diagonalized 



as H e ff = ^fC z where u s = a/ + d 2 e 2 . The transformed chirality operator is C z = 
U^C Z U = I fr)(fT | — | JJ-)(-D- | where the transformation operation is U = \ + \if)J)(ty \. 

The eigenvector of H eff are |^±) = ±^k\C = +1, M = \) + y/5^\C = —1, M = |). The 
coefficients 5± = (u 8 ± A so )/2u s . Then, the interaction Hamiltonian in this dressed state 
basis | ft~ (JJ-)) can be given by 

H'j(t) = A\t) ® B{t) + A{t) ®B\t). (8) 

Here A(t) = ^9^'^% and Bi (t) = 5 Q C Z + 6+e iu ' t C + - (Le-^O. where 6 = y/5^L. 
The expression of the time-convolutionless master equation in the dressed state basis is 
obtained as 

^ = -i[H eff + H L ,p(t)] + mt)) + NL\p{t)Y (9) 
Where H L = lm(T - T' )5 2 C 2 + £ 9=± Im(r ? - Y' q )5 2 q C\C q . The parameters Y q and T' q (q = 
0, ±) are determined by 



r„ = / dU *P \aA 2 ■ n j e^-"- quJs)it - tl) 



Jo j 

r; = I dt! V \gj\ 2 ■ (nj + i) e K— (io) 
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where the mean number for the jth mode of the thermal environment 

at T temperature. The Lindblad superoperator in Eq. (9) is given by 

mt)] = E im{t)[c m pci - \{cic m ,p}i (ii) 

m=z,i: 



where the decay rates are obtained as j z (t) = 2£QRe(r + r' ), 7+(t) = 2<5(j~Re(r + ) + 
2<5 ~Re(r / _) and 7_(t) = 25 c j"Re(r_) + 25 c }~Re(r / + ). The notation Im(Re) denotes imaginary 
(real) part of a complex parameter. The last term in Eq. (9) is very complicate, 



NL[p{t)\ = r • %8 + {C x pC. - pC_C z ) - 5 5_(C z pC + - pC + C z )} 
+ T + • [6 6+(C + pC g - pC z C + ) - 5 + 5_{C + pC + - pC + C + )] 

- T_ ■ [5 5-(C-pC z - pC z C.) + 5 + 5^(C-pC- - pC-CJ)\ 
+ r ■ i5 5+(C_pC z - pC z CJ) - 5 5_(C + pC z - pC z C + )} 
+ T' + ■ [6 S + (C z pC + - -pC+C z ) - 5 + 5-{C + pC+ - pC + C + )) 

- rl ■ [6 6- {C z pC- - pC-C z ) + 5+5- (C-pC- - pC-C-)) + h.c. (12) 

The notation h.c. represents Hermitian conjugate. The first part in Eq. (9) is the unitary one. 
Hi is the Lamb shift Hamiltonian and describes a small shift in the energy of the eigenvectors 
of H e ff. The Lamb shift Hamiltonian has no qualitative effect on the decoherence of the 



system and may be neglected 23J]. Meanwhile, according to the calculation of the effective 
spin-electric coupling, the characteristic time for the low-energy molecular magnet of Cu% is 
about t s = u^ 1 ~ 10ns [32| which is always much smaller than the correlation time of the 
thermal environment te ~ 1ms [5|. Under the condition of t s rg, the influence of the last 
;erm in Eq. (9) on the decoherence of the chiral states is very small and usually negligible 
l|. However, it is also noticed that the dynamics of the non-Lindblad operator NL[p(t)] is 
indispensable for the decoherence under the circumstance of r s ^ te- In the following, the 
parts of H e fj and Lindblad superoperator L[p(t)] are dominant in the second-order time- 
convolutionless master equation describing the dynamics of spin states in an electrically 
driven single molecular magnet. 



III. NON-MARKOVIAN DYNAMICS OF SPIN STATES 

An efficient way to describe the dynamics of the low-energy spin chiral states is to analyze 
the Bloch vector R{t) = (R x , R y , R z ) for the reduced density matrix p{t) of the single 
molecular magnet. The three components of the Bloch vector are defined as 



R j = T^lptyCj], (j = x, y, z). 
6 



(13) 



For any initial state .R(O) = (sin 9 cos 0, sin 9 sin 0, cos 9), the analytical expression of the 
Bloch vector at time t are obtained as 

R x (t) = e" r(t) • cos(uj s t + 0) sin 9 
R y {t) = e~ r{t) ■ sm(u s t + 0) sin# 

R z (t) = e" p W • {costf + fdh ■ e p ^[ 1+ (h) - 7 -(*i)]}, (14) 

Jo 

where r(t) = \ f*dt 1 [-y+(t 1 ) + 'y-(t 1 ) + 4'y z (t 1 )] and p(t) = /„* dfi[7+(*0 + 7-(*i)]- Th e angles 
satisfy that 6 1 G [0, 7r] and G [0, 2n]. The traditional Lorentzian spectral density function 
is used to describe the thermal bosonic environment such as the quantized electromagnetic 
field inside a cavity. The weak interactions between the system and thermal environment are 
given by the spectral density function, Ylj \dj\ 2 ~^ Jo°° Ji^^duj' where J(u/) = • 
The weak coupling constant u 2 <C u s has frequency dimension. The correlation time scale 
is obtained as r E = A -1 where A denotes the width of the distribution quantifying leakage 
of photons, ojq is the center frequency of the electromagnetic field in the cavity. 

The time evolution of the decay rates 7™(t), (m — ±, z) is plotted in Fig. 1 for different 
temperatures. The 7 m (t), (m = ±,z) is plotted in Figs. 1(a), 1(b) and 1(c) for small value 
of = 0.1. The values of 7 ± always oscillate between some positive and negative values 
with slightly damping rate. The amplitude of 7±(i) for T = is much smaller than that 
for T = 1. In non-Markovian quantum jumps formalism, negative values are regarded as 
the occurrence of reversed quantum jumps which can indicate the non-Markovian dynamics 
induced by the environmental memory. The memory effects describe the exchange of energy 
and information between the system and environment. The 'y m (t),(m = ±,z) is plotted 
in Figs. 1(d), 1(e) and 1(f) for large value of "^y^ = 10. The curve of T = is almost 
the same as that of T = I. That is, the influence of low temperatures on the dynamics of 
the decay rates is almost negligible. This is because that larger values of ^^f^ represent 
smaller effective coupling between the system and environment. Meanwhile, the effects of 
temperatures are also very weak. However, for a small value of wo ~ w as shown in Figs. 1(a), 
1(b) and 1(c), the oscillation of the decay rates at low temperature is dramatic in comparison 
with that of T = 0. It is seen that the impact of 7 2 on the decoherence is dominant in the 
the condition of r s <C te- For large value of U) °~ a; , the non-Markovian oscillation of 7^ occurs 
due to the negative values of j z . For long time limit, the values of the decay rate j z gradually 
approach to some steady values close to zero. This means that the non-Markovian dynamics 
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on the correlation time scale is reduced to the Markovian limit on the long time scale. 

The z-component of Bloch vector R z is plotted in Fig. 2 as a function of At. The phe- 
nomenon of the collapses and revivals of R z indicates that the rapid exchange of energy 
and information happens between the system and environment. It is shown that the de- 
coherence of spin states can be suppressed to a certain extent when the value of is 
large. To enhance the coherence of low-energy spin states, the frequency u of the external 
electric field and the width A of spectral density function can be decreased. Fig. 3 also 
demonstrates the effects of the electric field on the decoherence. For larger values of ^ M2 -, 
the decoherence can also be restrained. It is found out that the coherence can be improved 
by the manipulation of the frequency oj of the electric field when u is almost resonant with 
the transition frequency u so of the spin-orbit interaction. The increase of the strength of 
the electric field is also useful for the control of the decoherence. 

In respect to quantum information processing, the stability of the information stor- 
age needs to be analyzed when encoding qubits in single molecular magnet are electric- 
controllable in the thermal environment. From the perspective of von Neumann entropy, 
the pointer state 33], 34 1 can be defined as one initial state which becomes minimally entan- 
gled with the environment during the evolution. The study of the pointer state can help us 
to understand the effects of the decoherence on quantum information processing. The en- 
tropy for the reduced density matrix p(t) of the non-Markovian decoherence can be written 
as 

E{t) = -Tr[p(t) lnp(t)] = £ vdnv h (15) 

i=l,2 

where the eigenvalues of p are given by i>j = — ^— ^ — — . By means of Eq. (14), the entropy 
E(t) is easily calculated and dependent on the angle 9 of the initial spin state. 

For zero temperature of T = 0, the entropy E during the decoherence is plotted in Figs. 
4. It is seen that the values of the entropy are always increased with the time in Fig. 
4(a). The larger values of the entropy denote the more entanglement between the system 
and environment. Moreover, the time-dependent behavior of the entropy also shows the 
non-Markovian oscillation in Fig. 4(b). The revivals of the entropy are induced by the 
environmental memory effects. When the initial state is at 9 = n, the values of E almost 
remain the minimal ones in the evolution. Therefore, the pointer state for the thermal 
environment with T = is the state of 9 = 9 P = n which is almost the ground state 1^-)- 



S 



The effects of the environment on the pointer state at T = is very weak. In 



act, the 



selection of the pointer state is determined by the properties of the environment |33l . [34| . At 
low temperature of T = 1, the dynamics of the entropy is shown in Fig. 5. From Fig. 5(a), 
it is clearly seen that the values of the entropy for T = 1 are always increased more quickly 
than those of T = 0. This means that the single molecular magnet is easily entangled with 
the environment at T ^ 0. It is found out that the pointer state with the initial angle 
6 P are varied with the parameter For large value of u '°~ u ' > 10, the pointer state is 

approximately the ground state because of the very small effective coupling between the 
system and the environment. 



IV. DISCUSSION 



The decoherence of the low-energy spin states in an electrically driven single molecular 
magnet weakly coupled to a thermal environment is investigated. By means of the time- 
convolutionless non-Markovian master equation, the reduced density matrix for the spin 
states can be derived in the condition of t s <C t^. In regard to the Lorentzian environment, 
the oscillations of the decay rates between positive values and negative ones appear. This 
phenomenon indicates the memory effects from the non-Markovian environment. The rapid 
non-Markovian decoherence of the Bloch vector occurs due to the quick exchange of energy 
and information between the system and the environment. The decoherence can be efficiently 
suppressed by adjusting the electric field and the parameters of the environmental spectral 
density function. In quantum information processing, the selection of the pointer states can 
be determined by the properties of the environment. 
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Figure Captions 

Figs. 1 

The time evolution of the decay rates j m (t), (m = ±, z) is plotted as a function of the 
scaled time At when the parameters are ^ = 100 an< l ^f 2 = 0-4. For (a), (b) and (c), the 
parameter is '^^^ = 0.1 while for (d), (e) and (f), the parameter is = 10. The black 
solid line denotes the case of T = 1 and red solid one represents that of T = 0. 

Fig. 2 

The dynamics of the z component of the Bloch vector R z is plotted as a function of the 
scaled time At for T — 1, ^ = 100 and ^ = 0.4. The black solid line denotes the case of 
^""^ = 10 and the red one represents that of """^ = 0.1. 

Fig. 3 

The dependence of the decoherence of the spin states on the electric field is plotted for 
<£, _ ioq ; <£o__h _ q i anc j t = 1. The parameter ^ can be modified with the frequency u 
or the strength e of the electric field. 

Figs. 4 

The dynamics of the von Neumann entropy E is plotted at T = 0. (a) The initial state 
are changed with 9 and the parameters are = 100, =0.1 and = 0.9. (b) The E 
is plotted for 9 = (black line), 9 — | (red line) and 9 = ir (green line). The non-Markovian 
dynamics of the entropy is clearly presented. 

Figs. 5 

The dynamics of the von Neumann entropy at T = 1 is plotted for ^ = 100 and ^f 2 - = 0.9. 
(a) The initial state are changed with 9 when ^o-^ _ ^ pointer state represented 
by 9 P is plotted as a function of the environment parameter . 
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